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Abstract 

We prove that, under fairly general conditions, a properly rescaled 
determinantal random point field converges to a generalized Gaussian 
random process. 



1 Introduction and Formulation of Results 

Let E he a locally compact Hausdorff space satisfying the second axiom of 
countability, B — cr-algebra of Borel subsets and fi a cr-finite measure on 
{E, B), such that fi{K) < oo for any compact K G E. We denote by X the 
space of locally finite configurations of particles in E: X = = {xi)'^_^ : 
Xi E E Vz, and for any compact K C E •= i^i^i '■ Xi E K) < +00}. 

A o"-algebra of measurable subsets of X is generated by the cylinder sets 
= G X : = ''^j) where 5 is a Borel set with a compact closure 

and n G = {0, 1, 2, ... }. Let P be a probability measure on (X, JF). A 
triple {X,J^,P) is called a random point field (process) (see [DVJ], [Lel- 
Le3]). In this paper we will be interested in a special class of random point 
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fields called determinantal random point fields. It should be noted that 
most, if not all the important examples of determinantal point fields arise 
when E = Y['i=i^i (here we use the notation ]J for the disjoint union), 
Ei = R'^ or Z'^ and /X is either the Lebegue or the counting measure. We will, 
however, develop our results in the general setting ( our arguments will not 
require significant changes). 

Let dxi, i = 1, . . . ,n be disjoint infinitesimally small subsets around the 
Xj's. Suppose that a probability to find a particle in each dxi (with no 
restrictions outside of Y[i=i dxi) is proportional to HILi f^{dxi), i.e. 

P{4{dxi) = 1, i = 1, . . . , n) = pn{xu Xn)lJi{dxi) . . . n{dxn) (1) 

The function p„(a;i, . . . ,Xn) is then called the n-point correlation function. 
The equivalent definition is given by the equalities 



where Si, . . . , are disjoint Borel sets with compact closures, m > 1, rij > 

1, i = 1, . . . , m, ni H h = n. 

A random point field is called determinantal if 



where K{x, y) is a kernel of an integral operator K : L'^{E, dfi) — >• L'^{E, d/j,) 
and K{x, y) satisfies some natural regularity conditions discussed below. 
Such a kernel K{x, y) is called a correlation kernel. 

It follows from 2) and the non-negativity of the n-point correlation func- 
tions that K must have non-negative minors, and in particular if K is Her- 
mitian it must be a non-negative operator. In this paper we shall always 
restrict ourselves to the Hermitian case. 

Determinantal (also known as fermion) random point fields were intro- 
duced by Macchi in the early seventies (see [Mai], [Ma2], [DVJ]). A recent 
survey of the subject with applications to random matrix theory, statistical 
mechanics, quantum mechanics, probability theory and representation theory 
is given in [Sol]. Diaconis and Evans in [DEI] introduced a generalization of 
determinantial random point processes, called immanantal point processes. 

Let be a Hermitian, locally trace class, integral operator on L'^{E, dn). 
Suppose that we can choose a kernel K{x, y) in such a way that for any Borel 




Pn(xi, . . . , x„) = det(X(xi, X 



'i>)l<i,j<n' 



(2) 
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set B with compact closure 



Tt{KXB)^ I K{x,x)dx, (3) 
Jb 



where Xb denotes the multiphcation operator by the indicator of -B (=projector 
on the subspace of the functions supported in B). 
Since it is always true that 

Tr{KXB, . . . KXbJ - 

/ K(XI,X2)K(X2,X3) . . .K(Xn,Xi)dljL(xi) . . .dljL(Xn) (3') 

J BlX---xBn 

for n > 1 and Borel sets Bi, . . . ,Bn with compact closure, (3) implies that 
(3') holds for all n. 

(3) can always be achieved for E = M.'^ (see e.g. [Sol], Lemmas 1,2). 
Prom now on we will assume that both (2) and (3) are satisfied. 

The main goal of our paper is to study the behavior of linear statistics 

i 

for sufficiently "nice" test functions in a scaling limit. The moments of Sf 
can be calculated from (2). For instance, 

ESf^ J f{x)K{x,x)di,{x), (4) 



Var Sf^ J f{x)K{x,x)df,{x) - j f{x)f{y)\K{x,y)\'df,{x)dM. (5) 

Taking E — 'M} and K{x,y) — ^^^^^z^, a so-called sine kernel, we obtain 
a random point field well known in the theory of random matrices. It can 
be viewed as a limit n — > oo of the distribution of the appropriately scaled 
eigenvalues of n x n random Hermitian matrices with Gaussian entries (see 
e.g. [D], chapter 5). It was proven by Spohn in [Sp] (see also [So2]), that if K 
is the sine kernel and a test function / is sufficiently smooth and fast decaying 
at infinity, then Xli^-oo /(l^) ~^ I^oo f{^)dx converges in distribution to the 
normal law A^(0, |/(^)P " \k\dk), where / is the Fourier transform of /, 
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f{k) — J^^f{x)e '^'^'^^^dx. In other words we can say that the random signed 
measure 

oc 

Six - f ) - Ldx 

i=—oo 

converges as L — oo to the generalized self-similar Gaussian random process 
with the spectral density (see e.g., [Dobl], [Dob2] §3, [S], and, for the 
introduction to the theory of generahzed random processes, [GV]). The fact 
that the variance of the hnear statistics Y^°^_^ /(^) does not grow to infinity 
for Schwartz functions is the manifestation of the strong repulsiveness of the 
distribution of the eigenvalues of random matrices. Similar results for other 
ensembles of random matrices have been obtained in [DS], [Jol], [Jo2], [B], 
[BW], [So2], [W],[DE2]. The kernels appearing in these ensembles are, in 
some respect, very much like the sine kernel. In particular, the variance of the 
number of particles in an interval grows as a logarithm of the mathematical 
expectation of the number of particles. The following result was established 
by Gostin and Lebowitz for the sine kernel ([GL]): let / be an indicator of 
an interval, f = Xj, I = {a,b), then 

oo 

E J2 fi^i/L) = E{#{x, -.aLKx.K hL)) = L{b - a), 

i=—oo 

oo 

Var /(VL)) = -logL + 0(l), 

i=— oo 

and 

^(xi : aL < Xi < hL) — L{b — a) 
V^logL 

converges in distribution as L ^ cxd to the normal law A^(0, 1). The proof of 
the Gostin-Lebowitz theorem holds, quite remarkably, for arbitrary determi- 
nantal random point fields with Hermitian kernel. 

Theorem ([So3]) Let (X, JF, Pj^), L > 0, be a family of determinantal ran- 
dom point fields with Hermitian locally trace class kernels Kl and {/l}l>o be 
a family of Borel subsets of E with compact closure. Then ifYaiL {#{xi Xi e 
II)) oo, the normalized random variable ^('^'■'^'f^) converges in 

distribution to N{0, 1). 
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Here and below we denote by E^, Var^ the mathematical expectation 
and the variance with respect to P^. One can also establish a similar result 
for the step functions (finite linear combinations of indicators). 

Theorem Let {X, T . Pl) be a family of determinantal random point fields 
with Hermitian locally trace class kernels Kl and . . . , I^j'^^ } l>o be a 

family of Borel subsets of E, disjoint for any fixed L, with compact clo- 
sure. Then if for some ai, . . . , e 'M} , the variance of the linear statistics 
Z)i^-oo ^^^^ /i(^) = • ^tU) (^); grows to infinity in such a way 

that VarL(#jO)) = ©(Var^, (Yl^-oo fj-i^i))) f^'^ ^'^V 1 < J < ^; ^he Central 
Limit Theorem holds: 



Ei=i «f ^ ■ #/0) - El (Ei=i ■ #70)) 



Ar(0, 1). 



Reiiicirk 1 We use standard notations / = 0{g) and / = o{g) when ^ stays 
bounded or ^ — > 0. 

9 

Remcirk 2 The last theorem has been explicitly stated in [So3] only in the 
special case of the Airy and Bessel kernels and the kernels arising in the 
classical compact groups (see Theorems 1, 2, 4, 6), however the key Lemmas 
7 and 8 proven there allow rather straightforward generalization to the case 
of an arbitrary Hermitian kernel. A result close to our Theorem 6 from [So3] 
was also established by K. Wieand ([W]). 

We recall that a Hermitian kernel if (x, y) defines a determinantal random 
point field if and only if the integral operator K is non-negative and bounded 
from above by the identity, 

< if < Id (6) 

([Sol], Theorem 3). For the translation- invariant kernels K{x — y) and E — 
■^d is equivalent to < X(i) < 1, where 

The sine kernel K[x — y) = corresponds to K{t) = ^\ i ii(^)) 

the indicator of [— |, |]- It might be worth noting and actually is not very 
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difficult to see, that the logarithmic rate of the growth of Var {^{xi : < 
L)) for the sine kernel is the slowest among all translation- invariant kernels 
corresponding to projectors, K = Xb, for which inf(B) and sup{B) are the 
density points of B. For the generic translation-invariant kernel K{x — y) 
( K is not an indicator) Var (^(xj : \xi\ < L)) is proportional to Vol {xi : 
\xi\ < L) ~ E{4{xi : \xi\ < L)) ([Sol], section 3). 

In our main result we prove CLT for the linear statistics when the variance 
grows faster than some arbitrary small, but fixed, power of the mathematical 
expectation. 

Theorem 1 Let {X,J-',Pl), L > be a family of determinantal random 
point fields with Hermitian correlation kernels K^. Suppose that fL, L >0 
are bounded measurable functions with precompact support (i.e. supp{fi) has 
a compact closure for any L > 0), such that 



Var^ Sf^ — > oo 05 L — > oo (8) 

and 

sup\fL{x)\^o{YeiTLy, SL%u=0((Vari%)^), (9) 

for any e > and some 6 > 0. Then the normalized linear statistics 
-^=== converges in distribution to the standard normal law N{0, 1). 

As a very important special case of Theorem 1 one can consider fiix) :— 
f{TLx), where {Tl}, L G is a one-parameter family of measurable trans- 
formations Tl : E E such that Tj^^D has compact closure for any compact 
D. If for a sufficiently rich class of test functions / (e.g., continuous func- 
tions with compact support) (8), (9) are satisfied, and the rate of the growth 
of Var^ (Sfj^) is the same, 

Var^ (%)=i?(/)-Vi- (1 + 0(1)), 

where B{f) is some functional on a space of test functions. Theorem 1 implies 
that the random signed measure 

VP ( S{x-TLXi)-T4KL{x,x)dfx{x)U 

V=— oo / 
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converges as L oo to the generalized Gaussian process with the correlation 
functional B{f,f) = B{f) ( we denote by TL{KL{x,x)dfj,{x)) the image of 
the measure KL{x,x)dij{x) under Tl). 

Let us consider a Euclidean one-particle space E = 'R.'^ , a. one-parameter 
family of dilations 

Tl-.W^^ W^, Tlx = x/L, 

and a correlation kernel 

KL{x,y) = AL{x-y) + RL{x,y), (10) 

where 

\RL{x,y)\ < Q{Xabs+yabs): (H) 

Xabs = (kil, . . . , \xd\), Q e L2(R^) n L°°(R^) It follows from (6), (10) and 
(11) that < < Id, which implies < ^^(A;) < 1, < J^dALik) - 
{AL{k))'^dk = Al{0) - J^a \ AL{x)\'^dx —. al, and ctl = if and only if Al is 
an indicator. 

Theorem 2 Let the kernel Ki satisfy (10), (11) and there exist constants 
const, (T > and kl ^ oo as L ^ oo such that 

a as L — > oo, 

|^l(0)| < const, 

and 

\AL{x)\^dx 0. 



/ 

J\x\ 



' \x\>L/kl 

Then for any real-valued function f e L^(]R'^) nL^(M'^) the normalized linear 
statistics 

converges in distribution to the Gaussian random variable N (O, J^^ ifi^))'^ dx) . 

Remark 3 Theorem 2 says that under the stated conditions the random 
signed measure 



oo 



7^ 



6 - ^) - Al{0) ■ Ldx 



\i=—oo 
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converges to the white noise as L — * oo (for the definition of the white noise 
see e.g. [H]). Similar results hold in the discrete case. 

Let us now restrict our attention to the translation-invariant kernels 
K{x, y) = A{x — y). We will use the notation 

m(A) := j A{k) - A{k)A{k - \)dk. 

Observe that cr^ = m(0) and 

Var /(^^)) = / l/(A)rm(A)dA. (12) 

In particular 

Var = Vol ([-L, L]") ■ (m(0) + o(l)) . (13) 

It follows from (12) that the rate of the growth of the variance of Sf^ depends 
on the asymptotics of m(A) near the origin. In the next theorem we consider 
the degenerate case cr^ = in one dimension. 

Theorem 3 Let K{x, y) = A{x — y) be a translation-invariant kernel in M} 
and m(A) = |A|"(/9(A), where v^(A) is a slowly varying function at the origin 
and < q; < 1. Then for any Schwartz function f : ESf^ = ^^4(0) J f{x)dx, 
Var Sf^ = Li-"99(L-i) / \f{k)\'^\k\''dk{l + o{l)), and 



(Li-"(/?(L-i))i 



converges in distribution to N{0, f \f{k)\'^\k\°'dk. 

Remark 5 We recall that (/?(A) > is slowly varying at the origin if 
limA^o — 1 for any a 7^ (see [Se]). 

Remcirk 6 The result of Theorem 3 can be interpreted as the convergence 
in distribution of the random signed measure 

1 

to the self-similar (also called automodel in the Russian literature) general- 
ized Gaussian random process with the spectral density < a < 1. 
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Self-similarity means that the distribution of the process is invariant under 
the action of the renorm-group ^{x) ^{ax)a'^ , 7 = . The self-similar 
generalized Gaussian random process corresponding to « = is exactly the 
white noise (see Remark 3 above). It was proven by Dobrushin that the only 
self-similar random processes in M} are the ones with the spectral density 
|/c|", 0<q;<1. A self-similar generalized random process with the spectral 
density \k\ appeared in the Spohn's results ([Sp]) discussed above after the 
formulas (4), (5) (see also [Jol], [B], [So2]). For additional information on 
self-similar random processes we refer the reader to [Dobl], [Dob2],[S]. 

Example Let A be the indicator of U„>i [n,n -|- n~'^], (3 > 1. Then m(A) = 
, ,1-- 

const ■ |A| /^(l + 0(1)). On the other hand if the length /„ of the n— th 
interval [n,n + /„] decays sufficiently fast, say < In+i < /^"'"S e > than 
m(A) is not regularly varying at the origin. 

Finally we consider the case when A is the indicator of a union of 1 < 
^ < 00 disjoint intervals. It is straightforward to see that then m(A) = £\X\ 
near the origin. 

Theorem 4 Let A be the indicator of I, I = Ui^i[ai,bi], ai < 61 < 02 < 
b2 < ■ ■ ■ < ae < bg. Then for any Schwartz function f Xli^-oo / ( l") ^ 
^4(0) ■ L f{x)dx converges in distribution to N{0,£ ■ \f{k)\'^\k\dk) . 

The proofs of Theorems 1-3 will be given in the next three sections. The 
proof of Theorem 4 is the same, modulo trivial alterations, as the one given 
for the sine kernel in [So2]. 

The author would hke to thank P.Deift and Ya. Sinai for useful discus- 
sions, T. Shirai for making available the preprint [ST] prior to its publication 
and the referee for useful comments. 

2 Proof of Theorem 1 

We are going to prove Theorem 1 by the method of moments. Let us denote 
by Cn{Sf) the n*'* cumulant oi Sf. We remind the reader that for a random 
variable rj with all finite moments , the cumulants Cn{r]), n — 1,2, . . . are 
defined through the Taylor coefficients of the logarithm of the characteristic 
function : 

oo 

log E{e^pittr))) = ^C„(7?)(ii)7n!. 

n=l 
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We show that the n*'^ cumulant of the normahzed hnear statistics — % 

(Var Sf^)2 

converges to zero as L ^ oo for sufficiently large n {n > max(2(5, 2)). The 

Sf ~ESf 

Lemma 3 from the Appendix then asserts that all cumulants of ^ 

(Var Sf^)2 

converge to the cumulants of the standard normal distribution, which implies 
the weak convergence. 

We recall the lemma established in [So2] (see formula (2.7)). 

Lemma 1 

m — 1 {ni,...,n',-n)-ni^ ^nm='n.. (-^^j 

rti>l, i=l,...,m 

K{xi,X2)r''{x2)K{x2,X3) . . . P"'{Xm)K{Xm-,Xi)dll{Xi) . . . dll{x^) 



Using Lemma 1 we will be able to estimate the cumulants of Sf^. We 
claim the following result to be true. 

Lemma 2 Under the assumptions of Theorem 1 

CniSfJ = 0((Vari%)^+^), n>l, (15) 

where e is arbitrarily small. 

Proof of Lemma 2 It follows from (14) that Cn{Sf^) is a linear combination 
of 

J f2' ixi)KLixi, X2)fT{x2)KL{x2, X3) . . . fl"^ {Xm)K L{Xm, ^i) 

d^{x^) . . . d^{x^) = Tiifl^KLfl^KL . . . /r ^l), 
where rii>l, i = 1, . . . ,m, m > 1. 

We claim that each term is ©((Var^S'/^)^^'^). Indeed, if m = 1, then 
\Trf2KL\ = \Jf2ix)KLix,x)dfi{x)\ < \\h\[^^' J \h{x)\KL{x,x)d^{x) = 
IIMlSo-'^'^i^l, = 0((Varz.5;J^+^). 

If m > 1, represent Tr{f2^KLf2^KL ■ ..JI'^Kl) as a linear combination 
of TY{f2]LKLf±^^L^L ■ ■ ■ fZj^Ki), where we use the notations /+ = max(/, 0), 
/_ = max(— /, 0). Let us fix the choice of ± in each of the factors. Using the 
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1 1 

cyclicity of the trace and the inequahty \Ti{AB)\ < (Tr(^^*))2(Tr(SS*))2 
for the Hilbert-Schmidt operators ([RS], section VI. 6), we obtain 

mf±:LKLfT,LKL . . . IITlKl)] = \Tr{fZLKLfT,LKL . . . fl:LKLf^L)\ < 

ni "-1 ^2 1 ■"'2 ^1 

22. lil 1 



(16) 

The first factor at the r.h.s. of (16) is equal (again by the cychcity of the 

trace) to [Tr(/^^^/s:L/±'L^i^)] ^ (in particular we note that TT{giKg2K) > 
for non-negative gi, g2). 



Since 

cn2 \2 



+ ftLYK) - Tr((/-^ + fZ:L)K{fl]L + fZL)K) = 
Var (%i ) > 0, 



we have 



< TT{fl],K,fl^^,K,) < \ (Tr((/-^ + fT^.fK,) - Tri fl],K,fl],K,) 



1. 



T^(fT,LKLfT,LKL) j < 2l^((/^:^ + /rL)'^L)) = 0(TV(|MX^))o((Var^%)^) 
= 0((Var^%)^+^). 

(17) 

(n2 
{f^j^KLf^j^KL . . . 



(18) 

where C = K^fT^^KLf^^, D = U^^^K^^fl^Kr. . . . f^-^Kr^U^ Note 
that C > and Tr(C) = TT{f^\KLf^\K) = 0((Vari%)^+^) by argu- 
ments similar to (17). Using \Ti{CDD*)\ < Tr(C) • \\DD*\\ = Tr(C) • 
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([RS], Section VI.6) and < • where K = (EI^i "-i) - "-2, 

we conclude that (18) is 0((VarL Together with (16) and (17) this 
concludes the proof of the lemma. H 
Let us now apply Lemma 2 to estimate the cumulants of the normalized 
hnear statistics. We have 



1 = 0, 



and, for n > 2, 



It follows from the Lemma 2 and (19) that i ^t^^ 1 Soes to zero if 

V V J 

n > 25. 

Lemma 3 from the Appendix then implies that all cumulants of the nor- 
malized linear statistics converge to the cumulants of the standard normal 
random variable, and weak convergence of the distributions follows. 

Theorem 1 is proven. □ 



3 Proof of Theorem 2 

Let (E.dfi) be {M.^,dx) and Tlx — Consider a real- valued function / e 
L^(]R'^) n L'^{M.^). The mathematical expectations of 5"/^ is equal to 



ESf,= f f{x/L)KL{x,x)dx^ f f{x/L)AL{0)dx 

+ / f{x/L)RL{x,x)dx = AL{0)-L'^ [ f{x)dx+ [ f{x/L)RL{x,x)dx. 

jRd J J 
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By (11) the absolute value of the second integral is bounded by the sum of 
the integrals 

/ |/(±a;i/L, . . . , ±Xd/L) \ Q {2xu • • • , ^xa) dx 

1 1 
< /'(±xi/L, . . . , ±Xd/L)d^ Q\2x)d^ = O [hi) 

Therefore, 

ESf^ = Al(0) ■L'^- I f{x)dx + 0{L^). (20) 
The variance of 5*/^ is given by 

Var% = J f{x/L)KL{x,x)dx- J f{x/L)f{y/L)\KL{x,y)\''dxdy^ 

A^m^ j f{x)dx- J f{x/L)f{y/L)\AL{x - y)\'dxdy + r{L), 

(21) 

where 

r(L) = J f{x/L)RL{x,x)dx-2 J f{x/L)f{y/L)AL{x-y)RL{y,x)dxdy- 
J fix/L)f{y/L)\RL{x,y)\''dxdy = n{L) + r^iL) + rs{L) . 

It follows from the assumptions of the theorem that the second term at the 
r.h.s. of (21) is equal to 

L'^J \f{k)\'\X?{k/L)dk^L'\X\m J \f{k)\'dk-{l + o{l)) = 

L'^ J \AL{x)\^dx J f{x)dx{l + o{l)). 
Indeed, 

I J \fik)\'i\A^^{k/L) - \A^m)dk\ < 



'\k\>KL J\k\<KL 
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Since 

\\M^{t)\ = I J AL{k)AL{k-t)dk\ 
< J AL{k)dk = Ai(0) < canst 
we note that the first integral is bounded from above by 

const / \f{k)fdk — > as L — > oo. 

J\k\>{KLy/^ 

To deal with the second integral we estimate from above 

\\A^^ik/L)-\A^^{0)\< 

ALf{t){exp{27ritk/L) - l)dt\ 



1/ ^/ 

J\t\>L/KL J\t\<L/KL 

[ \AL\\t)dt + 0{l/,/lf^L)) 

J\t\>L/KL 

o(l) + 0(l/^/(«^))=o(l). 



Therefore 



Var % = (^^l(0) - J \AL{x)\'^dx^ L'^ J f{x)dx + o{L'^) + ? 
^2^d j f(^x)dx + o{L'^) + r{L). 
We claim that 

r{L)^o{L% 

Consider first ri(L). By (11) it is bounded by the integrals 
/ f{±xjL,...,±xa/L)Q{2x)dx 
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All of these integrals are estimated in the same way. For example, 
/ f{x/L)Q{2x)dx = L'^ [ f{x)Q{2Lx)dx = 

L'J f{x)Q{2Lx)A:^Q^,^,^^^ydx + L' J f{x)Q{2Lx)X^Q^,^,)^^y 
<L'\\Q\\ooJ f{x)X^Q^,j^,^^^ydx + L''-'^ J f\x)dx = o{L''), 
since 

l{x : Q{2Lx) > —=) — > 0. 
To estimate r3(L) we need to estimate the integrals of the form 

/ \f(x/L)\\f(y/L)\Q'(x + y)dxdy^L'' [ g{z/L)Q\z)dz, (24) 

where g{z) = / \f{x)\ \f{z - x)\XT^d{x)XT^d^{z - x)dx. 

Since g{z) is bounded, continuous, and zero at the origin, we have 

(24) = L^giQ) J Q\z)dz{l + o(l)) = 0(1"). 

Finally, 

|r2(L)| = I J f{x/L)f{y/L)AL{x - y)R{y,x)dxdy\ 

< J \f{x/L)\\f{y/L)\\A4x-y)\'dxdy ' 

1 

I \f{x/L)\ \f{y/L)\ \Rr^{y,x)\'dxdy ' = O (lI) o (lI) = o {L') . 

Combining the above estimates, we prove (23), which implies 

Var % = '^^^'^/ /^(^)^^(l + '^(l))- (25) 

If / is bounded, the Central Limit Theorem then follows from Theorem 1 
(compactness of the support of / is not needed since / e L^iMJ^) n L'^ 
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guarantees that all moments of Sjj^ are finite). The proof in the case of 
the unbounded / follows by a rather standard approximation argument. We 
choose > to be sufficiently large and consider a truncated function 

(fix), ii\f(x)\<N 

fix) = In, if/(x)>7v 

[-N iif{x)<-N. 

Observe that 

^i^Zl 71^) = aU" = 7' + 

can be made arbitrarily small by choosing N and L sufficiently large. 
Since 

~ N (o, [ f\x)dx 

aL2 V J\x\<N 

and 

lim / p[x)dx = I f'^{x)dx, 
'^-^°°J\x\<N J 

the result follows. 

Theorem 2 is proven. □ 



4 Proof of Theorem 3 

We now turn to the proof of Theorem 3. It is enough to establish that 

Var % = L^-''ifi{L-^) J \k\''dk{l + o(l)). (26) 

The result then will follow from Theorem 1. We have (see (12)) 
Var% = J \f{LX)\^L^m{X)dX = L J \f{k)\^m{kL-^)dk ^ 
L J |/>)r \krL-'^v{kL-^)dk = L^-XL-1) (27) 



(^(L-i) 
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It was proven by Karamata ([Kl], [K2]) that any slowly varying function at 
the origin can be represented in some interval (0, b] as 



where is a bounded measurable function on (0,6], such that ri{x) c 
as a; — > (|c| < oo), and e{x) is a continuous function on (0,6] such that 
e{x) ^ as X ^ 0. (for a modern day reference we refer the reader to [Se] , 
Theorem 1.2; of course a similar representation holds for (p also on some 
interval [b', 0) of the negative semi-axis). In particular 

uniformly in k on compact subsets of M} \ {0}, and the following estimates 
hold uniformly in k for sufficiently large L 

constiA;-" < ^{kL'^) / ^{L'^) < constsA;'", for 1 < A; < L, (30) 



constgA;"^ < (p{kL-^)/(p{L-^) < const4A;^ for < A; < 1, (31) 

where constj, i = 1, . . . 4, n > are some constants. 
The estimates (28)-(31) imply 

r i/»r \k\'^^^^dk r \krdk. 

Prom the other side, the integral over > L is o(l) since / is a Schwartz 
function and m is bounded. 

Theorem 3 is proven. □ 
Remark 7 We learned very recently that similar results to our Theorem 2 
have been independently obtained (in the discrete case) by Tomoyuki Shirai 
and Yoichiro Takahashi in the preprint [ST]. 



Appendix 

For the convenience of the reader we give here the proof of a rather standard 
fact. 
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Lemma 3 Let {qi} be a family of random variables such that Ci^tjl) = 
0, C2(?7i) = 1 and Cn{i]L) converges to zero as L oo for all n > N, 
where N < oo. Then linii^oo c„(?7l) = for all n > 2 and rjL converges in 
distribution to N{0, 1). 

1 

Proof Denote = max(|cj(?7i)|-? , 1 < j < — 1). It is clear that di > 1- 
Consider the random variable 

Since c„(?7l) = Cn{i]L)/dl we have \cn{fjL)\ < 1 for all n and Cnivi) 

for n > N. Consider {N — l)-dimensional vector (ci(?7l), . . . , CN-i{fiL))- Let 

(ci, C2, . . . , Civ-i) be a limit point. The Marcinkiewicz theorem (see e.g. [L]) 

states that if all but a finite number of cumulants of a random variable are 

non-zero then the random variable must either have a Gaussian distribution 

or be a constant. In both cases we have Cj — for j > 2. Therefore 
1 

dL — (c2(?7l))2 = 1 for sufficiently large L and Cn{r]L) — * for n > 2. 

L— >oo 

Convergence of the cumulants of 77^ to the cumulants of A'"(0, 1) is equivalent 
to the convergence of the moments which in turn implies convergence in 
distribution. 
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